Hamiltonian thermodynamics of (i-dimensional {d > 4) Reissner-Nordstrom anti-de 
Sitter black holes with spherical, planar, and hyperbolic topology 

GoiiQalo A. S. Dias* 
Centra Multidisciplinar de Astrofisica - CENTRA 
Departamento de Fisica, 
Instituto Superior Tecnico - 1ST, 
Universidade Tecnica de Lisboa - UTL, 
Avenida Rovisco Pais 1, 1049-001 Lisboa, Portugal 



Jose P. S. Lemos''' 
Centro Multidisciplinar de Astrofisica - CENTRA 
Departamento de Fisica, 
Instituto Superior Tecnico - 1ST, 
Universidade Tecnica de Lisboa - UTL, 
Avenida Rovisco Pais 1, 1049-001 Lisboa, Portugal 

The Hamiltonian thermodynamics formalism is applied to the general d-dimensional Reissner- 
Nordstrom-anti-de Sitter black hole with spherical, planar, and hyperbolic horizon topology. After 
writing its action and performing a Legendre transformation, surface terms are added in order to 
guarantee a well defined variational principle with which to obtain sensible equations of motion, and 
also to allow later on the thermodynamical analysis. Then a Kuchaf canonical transformation is 
done, which changes from the metric canonical coordinates to the physical parameters coordinates. 
Again a well defined variational principle is guaranteed through boundary terms. These terms 
influence the fall-off conditions of the variables and at the same time the form of the new Lagrange 
multipliers. Reduction to the true degrees of freedom is performed, which are the conserved mass 
and charge of the black hole. Upon quantization a Lorentzian partition function Z is written for the 
grand canonical ensemble, where the temperature T and the electric potential are fixed at infinity. 
After imposing Euclidean boundary conditions on the partition function, the respective effective 
action /* , and thus the thermodynamical partition function, is determined for any dimension d and 
topology k. This is a quite general action. Several previous results can be then condensed in our 
single general formula for the effective action /*. Phase transitions are studied for the spherical 
case, and it is shown that all the other topologies have no phase transitions. A parallel with the 
Bose-Einstein condensation can be established. Finally, the expected values of energy, charge, and 
entropy are determined for the black hole solution. 

PACS numbers: 04.60.Ds, 04.20.Fy, 04.60. Gw, 04.60.Kz, 04.70.Dy, 04.20.Ha, 04.50. Gh 

I. INTRODUCTION 

There are many approaches to calculate the entropy and the thermodynamics of a black hole. One can follow the 
original route where methods of field second quantization in a collapsing object are used to calculate the temperature T 
[H , and then uses the black hole laws to find the corresponding entropy f3| . Or one can use the Euclidean path integral 
approach to quantum gravity [3, [J] and its developments to obtain those thermodynamic properties [5, 
M, [U, E, [i3 (see also ^ and compare with @|). There are still other methods. The method we follow here is the 
one that builds a Lorentzian Hamiltonian classical theory of the gravity field and possibly other fields, and then obtain 
a Lorentzian time evolution operator in the Schrodinger picture. Afterward one performs a Wick rotation from real to 
imaginary time, in order to find a well defined partition function. The prescription implicit in this approach is, first, 
find the Hamiltonian of the system, second, calculate the time evolution between a final state and an initial state, i.e., 
between the final bra vector state < 32, S2, i2|, with metric (72 on a spatial hypersurface S2 at some generic prescribed 
time t2, and the initial ket vector state, >, with metric gi on a spatial hypersurface Si at some generic 

prescribed time ti, and, third, EucHdeanize time. Here, the amplitude to propagate to a configuration < g2,'S2,t2\ 
from a configuration \gi,Tji,ti > , is represented by < 52, S2, ^2! exp (— ii/(<2 — ii)) |<7i, Si, ti > in the Schrodinger 
picture. Euclideanizing time, t2 — ti = —if3 and summing over a complete orthonormal basis of configurations gn one 
obtains in general the partition function Z = ^ exp [—(} {En — (Zn^)], of the field 5 at a temperature T = and at 
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some chemical potential 4>, where £"„ is the eigenenergy of the eigenstate gn, and g„ is the eigenvalue of some variable 
conjugate to (j), such as the particle number (when (j) is effectively the chemical potential), or the charge number 
(when (j) is the electric potential), for instance. This route is based on the Hamiltonian methods for covariant gravity 
theories (isl. fl^. fl7| and for black hole spacetimes on the Hamiltonian method of Kuchaf It was developed by 
Louko and Whiting in (l9l | for the specific problem of finding black hole entropies and thermodynamic properties. 

Indeed Kuchar applied his procedure to the full vacuum Schwarzschild black hole spacetime, which is really 
composed of a spherically symmetric white hole plus a black hole plus two asymptotically flat regions. These regions 
are well pictured in a Carter-Penrose diagram. By considering the spacelike foliations of the full manifold, the true 
dynamical degree of freedom of the phase space of the Schwarzschild black hole was found. This degree of freedom, 
represented by one pair of canonical variables, is composed of the mass M of the solution and its conjugate momentum, 
which ph ysic ally represents the difference between the Killing times at right and left spatial inflnities. Louko and 
Whiting [19(1, by adapting the formalism to a spacelike foliation to the right of the future event horizon of the 
solution, and imposing appropriate boundary conditions, obtained the Hamiltonian H , the time evolution operator 
exp{—iHt) in the Schrodinger picture, the associated partition function and flnally the thermodynamics. This method 
has been applied for various theories of gravity, in several different dimensions, and in either asymptotically flat or 
asymptotically AdS spacetimes. 

There are several applications of the formalism to four-dimensional spacetimes (l9l. [20l. [2H . Now, it is clear that it is 
important to understand the physics in several different dimensions. Indeed, hints from many places like string theory, 
AdS/CFT (Anti de Sitter/conformal fleld theory) conjecture, extra large dimensions and the connected braneworld 
scenarios, point out to the possibility of a world with other space dimensions. A first attempt to study the thermo- 
dynamics of a higher dimensional black hole using the Hamiltonian approach was done for d = 5 Lovelock gravity 
[2a |. A related formaHsm developed in jH, US] studied some higher dimensional black holes that admit a reduction 
to a two-dimensional dilaton-gravity theory. There are also incursions in the app lication of the formalism into lower 
dimensional dilaton-gravity theories, e.g., in d = 2 [2BI, and in d = 3 [2^ . l27l |. In [2^ a connection between the 
Hamiltonian formalism of [3| and the path-integral formalism of @ is made. 

So here, we are interested in applying the Hamiltonian formalism introduced by Louko and Whiting (see 
also [23,[2i|,[23,[23,[2i,[2l,[2i,[23,[23|), to the d-dimensional, d > 4, Reissner-Nordstrom black holes with spherical, 
planar, and hyperbolic horizon compact topology, and with a negative cosmological constant, i.e., in an asymptotically 
anti-de Sitter (AdS) spacetime. Spherically symmetric electrically charged black holes in d dimensions are also known 
as Tangherlini black holes [29,], or simply d-dimensional Reissner-Nordstrom AdS black holes, see also [SQ] for the 
Kerr the d-dimensional Kerr black holes. Planar toroidal compact black holes in four dimensions were discussed in 
[sil. [s^l and charged ones in [s^l (see also [s^. [35|l. Hyperbolic toroidal compact black holes in four dimensions were 
studied in jsH] , and the charged version of the three together in [2l| . In d dimensions these topological black holes 
were analyzed in see also [H, [3§|. The solution for Kerr- AdS black holes in d dimensions was found in [iOl- One 
notes that if the cosmological constant is nonnegative, black holes with nonspherical topology do not exist. 

To develop the Hamiltonian formalism for these d-dimensional AdS black holes an important issue one has to 
deal with is to find the asymptotic fall-off conditions at infinity. In four dimensions the problem was settled by 
Henneaux and Teitelboim [41j, where through the precise asymptotic structure of the Kerr- AdS metric and acting 
on that structure with the four-dimensional AdS group one was able to find the correct fall-off conditions. Louko 
and Winters-Hilt [l^l used these fall-off conditions to study the four-dimensional Reissner-Nordstrom- AdS black hole. 
In d dimensions no such procedure is available, although the Kerr-AdS metric solution in d dimensions has already 
been found Thus, by following the same procedure as in [4l| one should be able to arrive at the correct fall-off 
conditions. Fortunately, in our case the problem is much simpler and we can adapt with some care to d dimensions the 
fall-off conditions of [2y|. Then, after following other procedures, one finds the reduced Hamiltonian, and a statistical 
analysis can finally be performed. We study the grand canonical ensemble and find a very general effective action 
/*, or equivalently, a very general partition function. From the action /* one can readily extract all the relevant 
thermodynamical information. We study the phase transitions of the system and the thermodynamic quantities, 
such as temperature and entropy, of the solution containing a black hole. The phase transitions are between hot fiat 
space and black holes, two different sectors of the solution space, one with trivial topology, the other with black hole 
topology, respectively. These phase transitions show similarities with the Bose-Einstein condensation phenomenon. 

There are several particular results that can be obtained from our general action formula /*. First, within works 
using a Hamiltonian thermodynamics formalism we recover the thermodynamics for spherically symmetric Reissner- 
Nordstrom- AdS black holes in four dimensions {d = 4) the thermodynamics for AdS black holes with planar and 
hyperbolic topology in four dimensions {d = 4) [2l|, and the results for the spherical black hole in a finite box with 
radius tb (when tb 00) in a Lovelock d = 5 theory, which can also be called a Gauss-Bonnet theory, without a 
cosmological constant Second, within works using an Euclidean path integral method we also recover from our 
general action formula, when appropriate, the thermodynamics of the d-dimensional spherical Reissner-Nordstrom 
black hole studied in (l2 . [l3| , as well as the thermodynamics of the d-dimensional Schwarzschild black hole studied in 
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[ll|. Moreover, we obtain the thermodynamics of the Reissner-Nordstrom-AdS black holes in four dimensions (d = 4) 
with spherically symmetric studied in which in turn recovers results found in 0,0,0] and in as well as the 
thermodynamics obtained in for the charged AdS black holes in planar topology. 

The paper layout is thus the following: in Sec. [ll] we write the ansatz for the black hole solutions of the charged 
black holes in d dimensions with spherical, planar toroidal, and hyperbolic toroidal topology, plus the ansatz for the 
vector potential A, following it with the metric and vector potential solutions for the d-dimensional black holes. The 
Carter-Penrose diagram of the charged black holes in d dimensions with spherical, planar toroidal, and hyperbolic 
toroidal topology, is then depicted. In Sec. the ADM form for the metric and A are spelled out, the canonical 
description for the Einstein-Maxwell action in d spacetime dimensions is done and is followed by the reconstruction of 
the action. After the canonical transformations, we reduce the Hamiltonian to the physical degrees of freedom. There 
follows the canonical quantization through the Schrodinger representation of the time evolution operator. In Sec. 
IIVI the thermodynamics is studied. First, there is the construction of the partition function for the grand canonical 
ensemble through the imposition of Euclidean boundary conditions on the time evolution operator found previously. 
Then, the critical points of the effective action in the partition function integral are obtained. Next, the effective 
action is evaluated at the critical points, where, depending on the value of the effective action at these critical points, 
one can obtain hot flat space or a black hole through a phase transition. The results are compared to previous results. 
It is then possible to determine the expected values of the energy and charge, plus the entropy, from the saddle point 
approximation of the partition function of the ensemble. In Sec. |V]we conclude. We choose G — 1, c — 1, h ~ 1, and 
/cb = 1 throughout. 



II. CHARGED BLACK HOLES IN d DIMENSIONS WITH SPHERICAL, PLANAR, AND 
HYPERBOLIC TOPOLOGY IN ASYMPTOTICALLY ADS SPACETIMES 

A. Solutions with spherical, planar, and hyperbolic compact topology 

The electrically charged AdS black hole solution can be generically written as (see [39]) 

ds^ = -F{R) dT^ + F{R)-^ dR^ + R^{dn^^_2f, (1) 

where we have d > 4, T is the Schwarzschild time coordinate, R is the Schwarzschild time coordinate, and F{R) 
is a function that yields horizons and thus black holes. Its form depends on the theory, we will be interested in 
d-dimensional general relativity. The constant k has values k = 1,0,-1 whether the topology is spherical, planar 
toroidal, or hyperbolic toroidal, all three being compact topologies. The angular part of the metric |(T]) is, for each 
k = 1,0,-1, 

d-3 

{d^i\_2f^del + sin^ eidel + ---^Y\ sit? e, del_2 , 
{<m\_2f=dB\ + ddl + del + --- + ddl_, , 

d-3 

{dil-^^f^dej + sinh^ 0idel+--- +sinh^ 61]^ sin^ 9, del_2 . (2) 

1=2 

In the spherical case, fc = 1, one can take the range of coordinates as the usual one < 9i < tt, < 9i < 2 n for i > 2. 
In the planar case, fc = 0, the range is arbitrary, though finite for a compact surface, and we can choose < 6*^ < 2 tt. 
In the hyperbolic-compact surface case, fc = — 1, the situation is in general more involved, and we need to choose a 
finite area surface for the black hole horizon in order to be able to study its thermodynamics. For example, in d = 4 
spacetime dimensions, the horizon is a 2-dimensional hyperbolic torus, i.e, a Riemann surface with genus (/onus > 2. 
For more details, see [3l| and [s^] (see also [s^). The vector potential one-form which is generically written as 

A = A{R) dT , (3) 

where the electric potential A{R) depends on the theory one is studying. We will work with general relativity in 
d-dimensions. 
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B. Black hole solutions: metric and vector potential 



We are interested in general relativity coupled to a Maxwell field in d-dimensions whose action is 

S^J d'xV^ [R+{d- l)id - 2) - ^F^.F^^'^^ , (4) 

with g being the determinant of the metric g^^, R is the Ricci scalar constructed from the Riemann tensor RafSjS, 
I is the AdS length, which the makes up the cosmological constant through — 3/~^, and Ff^^ = d^Ay — d^A^ is the 
Maxwell tensor, with A = A^dx^ the one-form vector potential. 

The electrically charged AdS black holes solutions of the Einstein-Maxwell- AdS action in d dimensions, have F{R) 
given by 

ni?) = fc + r^i?^-|^ + ^, (5) 

Here, the quantities M and Q are the mass and the charge parameters, respectively. The definition of mass in general 
relativity is always ambiguous. Usually one defines the Arnowitt-Deser-Misner (ADM) mass m, which in principle, for 
asymptotically well defined spacetimes is a well defined quantity. The black holes we study here are asymptotically 
AdS and so have a well defined mass. However, depending on the values one gives for the cycHc coordinates of the 
torus, one gets different answers in the planar and hyperbolic toroidal cases, not so in the spherical case which has 
well defined angular coordinates. Moreover, when performing our canonical analysis, as it will be done here, a mass 
which is the one necessary to make the important canonical transformations, and renders the expressions most simple, 
pops up. We call this mass the canonical mass M. The relation between the canonical mass and the mass parameter 

is, as we will see, 

= 8^ Af, where S^_2 is the area of the [d — 2) unit surface. For instance, in the spherical 
case and for zero cosmological constant, one has M = (d-2) s ^ ' where now in the spherical case M is also the ADM 

mass, M = and S/_2 is the area of the {d—2) unit sphere Q|. The same applies for the electrical charge, although 
here the ambiguities are not as strong. Here the canonical charge Q, i.e., the one necessary to make the important 
canonical transformations, and the charge parameter are equal, Q = Q. None of these is the ADM (Gauss) charge 

q. For instance, in the spherical case and for zero cosmological constant, one has — (rf-2)|d-3) ) • '^^^ vector 
potential function A{r) is 



where the relevant quantities have been defined above. 

The corresponding Carter-Penrose diagram is given in Fig. [TJ 




FIG. 1: The Carter-Penrose diagram for the Reissner-Nordstrom-AdS black hole with spherical, planar toroidal, and hyperbolic 
toroidal topologies, where is the outer horizon radius, Ri is the inner horizon radius, and the double line is the timelike 
singularity at -R = 0. The static region I, from the bifurcation (d — 2)-manifold to spacelike infinity R — oo, is the relevant region 
for our analysis. The regions II and 11' are inside the outer horizon, beyond the foliation domain. Region I' is the symmetric of 
region I, toward left infinity. Note that any point in the diagram is a (d — 2)-dimensional manifold. The bifurcation manifold 
can be a sphere, a planar torus, or a hyperboloidal torus, always (d — 2)-dimensional. 
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III. HAMILTONIAN THERMODYNAMICS FORMALISM 
A. ADM form of the metric and vector potential 

The ansatz for the metric field with which we start our canonical analysis is given by 

rfs2 ^ ^)2^^2 ^ _^ ^^^^^2 _^ ^(^^ ^^2 f^d-2)^ > (7) 

where t and r are the time and radial ADM coordinates used in the ADM metric ansatz for the black hole solutions 
of the Einstein-Maxwell-AdS action. In the subsequent developments we follow the basic formalism developed by 
Kuchaf The canonical coordinates are R and A which are functions of t and r, R = R{t,r), A = A{t,r). Now, 
r = is generically on the horizon as analyzed in [18|, but for our purposes r = represents the horizon bifurcation 
{d - 2)-manifold of the Carter-Penrose diagram (see, [H, [10, HU, iH, El, Hi, HI, S, [13, Ell). The coordinate r 
tends to oo as the coordinates themselves tend to infinity, and t is another time coordinate. The remaining functions 
are the lapse N — N{t,r) and shift functions — N^{t,r) and will play the role of Lagrange multipliers of the 
Hamiltonian of the theory. The canonical coordinates R — R{t,r), A = A{t,r) and the lapse function N = N{t,r) 
are taken to be positive. The angular coordinates are left untouched, due to the symmetries assumed. The ansatz 
([7]) is written in order to perform the foliation of spacetime into spacelike hypersurfaces, and thus separates the 
spatial part of the spacetime from the temporal part. Indeed, the canonical analysis requires the explicit separation 
of the time coordinate from the other space coordinates, and so in all expressions time is treated separately from 
the other coordinates. It breaks explicit covariance of the Einstein-Maxwell-AdS theory, but it is necessary in order 
to perform the Hamiltonian analysis. The metric coefficients of the induced metric on the hypersurfaces become 
the canonical variables, and the momenta are determined in the usual way, by replacing the time derivatives of the 
canonical variables, the velocities. Then, using the Hamiltonian one builds a time evolution operator to construct an 
appropriate thermodynamic ensemble for the geometries of a quantum theory of gravity. 
The canonical description of the vector potential one-form is given by 

A = Tdr + ^dt, (8) 

where F and $ are functions of t and r, i. e., r(t, r) and <i>(t, r). The function r(t, r) is the canonical coordinate 
associated with the electric field, and the function ^{t, r) is the Lagrange multiplier associated with the electromagnetic 
constraint, which is Gauss' Law. 



B. Canonical formalism 



We now replace the ansatz for the metric ^ and the ansatz for the one-form vector potential ^ into the action 
in Eq. (0]), obtaining 

S^[A, R,r;N, N'', ^] = J dt J dr a \k B N AR'^-'^ + Qr'^N AR'^~'^ ~ 2{d - 2)N-^ RAR'^~^ 

-BN-^AR^R'^-^ + 2{d - 2)N-^R {AN'')' + 2{d - 2)N-^A {R'N'') 

-2{d - 2)N-^ {AN'-y {R'N'') R"^-^ + 2BN-^ R' RN'' AR'^-'^ - BN'^A {R'N^'f R'^'^ 

-2{d - 2)N {A-^y R' R"^-^ - BNA-^ {R'f R'^-^ - 2{d - 2)NA-^R"R'^-^ 

+8TrN-^A-^ (f - $') ^ R'^-^ I , (9) 



where B is defined as B — (d ~ 3){d ~ 2), a is defined as 



and, e. g., for k — 1 we have S^_2 = Stt^'' ^)/^/r(^^), r(a::) being the Gamma function. 
The conjugate momenta are 



Pa = -2a{d-2)R'^-^N-^ (^R- N^'R'j , (11) 
Pr = -2a{d- 2)R'^-^N-^ \r (a - (A7V)') +{d- 3)A (r ~ TVi?')] , (12) 
Pr = IQ-naN-^A-^ {t -^'^ R"^-^ . (13) 
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After a Legendre transformation we find 



Sy 



A, R, r, Pa, Pb. Pr;N, N^, $ 



dt 



dr PaA + PnR + PrT - NH - N''Hr - 



•) , (14) 



where <& is defined as $ = <f> — N'T, and the constraints are 

(d — 3) . „T 1 



H = 



^p2u-(d-2) 



2 aid -2) 



Aa{d-2) IGira 

+ a (-kBAR"^-^ - Ql-'^KR'^-^ + 2{d - 2) (A-^)' R'R 



PrPaR 



+ BNA-^ {R'y R'^-^ + 2{d - 2)NA-^R"R'^-^ 

Hr = PrR' - AP't^ - VP^ , 
G = -P^. 



(15) 

(16) 
(17) 



The equations of motion are obtained through the variation of the action (fT4l) . with the constraints defined in Eqs. 
(JUl-dlZl), i.e., 



A 
R 

r 

Pa 



Pr - 



(2 a{d - 2))-^NAPr^R-^'^~^^ ~ (2 a{d - 2))-^ N PrR-^'^-^^ + (TV A)' , 
-(2 a{d - 2))-^NPkR-^'^-^^ + N^'R' , 
{Snay^NAPrR-'-'^-^'^ + (NT)' + $' , 

^l-i-Arp2^-('i-2) _ J_NP^R-^^-^^ +akNBR''-^ 

4a{d-2) ^ 167ra ^ 

+6r 2 a NR'^''^ - [2 a{d - 2)NR'R'^^^] ' A^^ a NB (R'f A'^R 

+2 a{d - 2)NR"A-^R'^-^ + N^'P^ , 

3 



2 



(18) 
(19) 
(20) 



(21) 



''^-NAPlR-^"-^^ + ^J^NAP^R-^"-^ 



-NPrPaR 



-(d-2) 



Pr 



4a '"'^ ' IGna 2a(d-2) 

+ak{d~ 4) BNAR'^-^ + 6r^ a{d - 2)NAR'^-^ +\2aBN (A^^)' 

-a{d- 4)BNA-\R')^R'^-^' + [2aBNA-^R'^-^]' R' 
- (2 a{d - 2)NA-^R'^-^)" + {N"" Pr)' , 



R 



d-3 



(22) 
(23) 



In order to have a well defined variational principle, with which to derive the equations of motion above (fT8l) - ((23l) . we 
need to eliminate the surface terms which remain from the variation of the original bulk action. These surface terms 
are eliminated through judicious choice of extra surface terms which should be added to the action lfT4|) . The action 
(fT4l) has the following extra surface terms, after variation 



Surface terms = 2 aNA'^ {R'^-^) SA + 2 a{d ~ 2)N' A-^ R'^-^SR -2a{d- 2)NA-'^R'^-^6R' 

-N^'PrSR + N'^ASPa + NTSPr + ^SPr . (24) 
In order to evaluate this expression, we need to know the asymptotic conditions of the above functions individually. 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 
(31) 
(32) 
(33) 



A(t, r) 


= Ao(t) + 




R{t, r) 


= Po(t)4 


R2{ty + 0{r^), 


PA{t, r) 


- 0(r3), 




PR{t,r) 


= 0(r), 




N{t,r) 


= Niity 










r(t,r) 


= 0{r), 




Pr{t,r) 


= Xo{t)^ 




$(i,r) 


= + 
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It is useful to red efine Xp and X2 in terms of two new quantities Qq and Q2 , such that Xq = Kd Qo and X2 = Kd Q2 , 
with Kd = Sna^y {d — 2){d — 3)/2n. The surface terms for r ^ suffer a modification in the definition Pr, in relation 
to the d = 4 case (2ol |. 

For r ^ 00 the fall-off conditions have to deal with the AdS asymptotic properties in d dimensions. In d = 4, the 
results in j4l| were the basis for the r ^ 00 fall-off adopted in [20|. Since one now knows the Kerr- AdS solution in 
d dimensions ji^l one could follow [4l!| to find the fall-off conditions for our case. Here we take a more pragmatic 
approach, by simply generaHzing, with some care, to d dimensions the results in |2fl| for c? = 4. Hence, we assume the 
following conditions: 

-1)), (34) 

(35) 
(36) 
(37) 

(38) 

(39) 
(40) 
(41) 
(42) 



A(t, r) 


= filr-') + fXit)r- 


d + 0°°ir 


R{t, r) 


= r + fpity-^"--"^ 4 






= 0°^ir- 






PB.it, r) 


= 0°^ir- 






N{t,r) 


- i?'A-i 


N+{t)+0 


oO('j.-(d+l) 


iV(t,r) 


= 0°^ir- 






nt,r) 


= 0°°ir- 






Prit,r) 


= x+it)- 






^it,r) 


= *+(^)^ 


-O^ir-^'^- 





where we defined the function f{lr ^) as 

/(/r-i) -Ir- 
with the double factorial being defined by 



E 



(2s- 1)!! 



2s 




1 s > odd 

2 s > even 
s = -1,0 



(43) 



(44) 



where [x] is the integer part of a given x E Z/2. Again, / is the AdS length, and k yields the topology of the horizon. 
As in r ^ case, it is also useful to define X+ = Kd Q+, with Kd defined above. 

The surfaces terms which must be added to ifTj)) in order to obtain a well defined variational principle are then 



Sdi 



A, R, Xo, X+;N, $0, 3'^ 



dt (2 aNiA^^ R"^-^ 



N+M+ 



where M+(i) is defined as 



M+{t) = 2 aid- 2) (A(t) -id - l)p(t)) 



(45) 



(46) 



One also should impose the condition of fixing A^iAq ^ and $0 on the horizon r — > 0, which is the same as saying 
that 5(7ViAq ^) = and — 0. The first condition fixes the rate of the boost suffered by the future unit normal 
to the constant t hypersurfaces at the bifurcation horizon and is important to have a well defined Hamiltonian 
thermodynamics. The second condition fixes the electric potential at the horizon. At infinity, r ^ 00, one fixes N+ 
and $+, which means 5N+ = and 5$+ = 0, in order that the variational principle is completely defined. So the 
total action is then the sum of the integral in Eq. (|45l) and the integral Eq. lfT4|) . 



E+as 



A, i?, r, Pa, Pr., Pr;N, iV^ $ 



Sy 



A, R, r. Pa, Pr, Pr;N, N^, $ 



A, R, Xo, X+;N, $0, *h 



(47) 



With this action one then derives the equations of motion (|T8|1 - ([23| without extra surface terms, provided that the 
fixing of 7ViAo\ ^0, N+, and $+ is taken into account [1^, [H, [13, Si, HI, S [13, [111 . 
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C. Reconstruction, canonical transformation, and action 

In order to reconstruct the mass and the time from the canonical data, which amounts to making a canonical 
transformation, we have to rewrite the form of the solutions of Eqs. ([l]). In the process, we have to consider how to 
reconstruct the charge from the canonical data, which is on the hypersurface embedded in this charged spacetime. We 
follow Kuchaf for this reconstruction. We concentrate our analysis on the right static region of the Carter-Penrose 
diagram. 

Developing the Killing time T as function of {t, r) in the expression for the vector potential |[6]), and making use of 
the gauge freedom that allows us to write 



where ^(i, r) is an arbitrary continuous function of t and r, we can write the one-form potential as 



Equating the expression l(49l) with Eq. ^ , and making use of the definition of Pr in Eq. (fT3)) , we arrive at 

Pr = KdQ- (50) 



with Kd = Aa^2i: (d — 2)(d — 3). We also obtain the derivative of the gauge function with respect to r 

^' = r + R-^R-^'-'^-'^^PrF-^APA . (51) 
In the static region, we have defined F as 

F{t,r)^k+l-^R^-^ + j^. (52) 
We now make the following substitutions 

T^T{t,r), R = R{t,r), (53) 

into the solution |(T|), getting 

ds^ = -{FT^ -F-^R^)dt^ + 2{-FT'T + F-'^R'R)dtdr + {-F{T'f + F~'^R^)dr'^ + R^{dn^d~-2f ■ (54) 

This introduces the ADM foliation directly into the solutions. Comparing it with the ADM metric |[7]), written in 
another form as 

ds^ = -(iV2 - K^{N''f) dt^ + 2K^N'' dtdr + K^dr"^ + R^ {dnd-2f , (55) 
we can write a set of three equations 

A2 = -F{T')'^ + F-\R')\ (56) 

A^iV^ -FT'T + P-^R'R, (57) 

_ A2(iV'^)2 ^ Ff^-F^^R^. (58) 

The first two equations, Eqs. l(56|) and Eq. l(Ff|) . give 

^ -FTt + F-^R'R 

-F{T'Y + F-^{R'f ^ ' 

This one solution, together with Eq. (|56|) . give 

R't-T'R , , 

N ^ , . 60 

^-F(r')2 + F-i(i?')^ 
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One can show that N(t,r) is positive (see Next, putting Eqs. l(59|) - (|60|) . into the definition of the conjugate 

momentum of the canonical coordinate A, given in Eq. (flTI) . one finds the spatial derivative of T{t, r) as a function of 
the canonical coordinates, i.e., 

- T' = (2 a{d - 2))-^R-^'^-^^F-^APa . (61) 

Later we will see that —T' ~ Pm, as it will be conjugate to a new canonical coordinate M, defined below in Eq. 
Following this procedure to the end, we may then find the form of the new coordinate M{t, r), as a function of t and 
r. First, we need to know the form of F as a function of the canonical pair A , R. For that, we replace back into Eq. 
((56)) the definition, in Eq. ((6T|) . of T' , giving 



A J \2a{d-2)Rd 
Equating this form of F with Eq. (|52|) , we obtain 



1 / P'^K^^ 



M{t, r) = -R'^-' \k + l-'R' + -^^^ -Fj , (63) 

where F is given in Eq. l(62|). However, it will be more convenient for what is coming to define the following mass 

M{t,r) = 2a{d-2)M(t,r) , (64) 

where M(t, r) is defined in Eq. ((63|) . The new canonical coordinate is thus M, which can also be called the canonical 
mass. It is now a straightforward calculation to determine the Poisson bracket of this variable with Pm = —T' and 
see that they are conjugate, thus making Eq. l(6T|) the conjugate momentum of M, i.e., 

Pm = (2 a{d - 2))-^ R-'^'^-^^ KPj^ . (65) 

In the same fashion, one needs to find the charge canonical coordinate related to Q. In this case the relation is trivial, 
and one finds that the canonical charge coordinate Q, or simply the canonical charge, is given by 

Q = Q. (66) 

Then, the other natural transformation is the one that relates Pr with the canonical charge Q. Through the relation 
of Pr with the charge parameter Q, Pr = KdQ, the canonical transformation between the old variable Pr and the 
new canonical variable Q is then immediately found to be Pr = KdQ. Like M previously, Q is a natural choice for the 
canonical coordinate, because it is equal to the charge parameter of the metric solution in Eq. H]). The one coordinate 
which remains to be found is Pq, the conjugate momentum to the charge. It is also necessary to find out the other 
new canonical variable which commutes with M, Pm, and Q, and which guarantees, with its conjugate momentum, 
that the transformation from A, R, F, to M, Q, and the new variable is canonical. Immediately is it seen that R 
commutes with M, Pm, and Q. It is then a candidate. It remains to be seen whether Pr also commutes with M, Pm, 
and Q. As with P, it is straightforward to see that Pr does not commute with M and Pm, as these contain powers of 
R in their definitions, and {P(<, r), PR(t,r*)} = 5{r — r*). So, rename the canonical variable P as P = R. We have 
then to find a new conjugate momentum to R which also commutes with M, Pm, and Q, making the transformation 
from {A, P, F; Pa, Pr, Pr } ^ {M, R, Q\ Pm, Ppj^, Pq } a canonical one. The way to proceed is to look at the 
constraint i/^, which is called in this formalism the super- momentum. This is the constraint which generates spatial 
diffeomorphisms in all variables. Its form, in the initial canonical coordinates, is H,. = —A P^ + Pr R' — TPy- In this 
formulation, A is a spatial density, P is a spatial scalar, and F is also a spatial density. As the new variables, M , 
R, and Q, are spatial scalars, the generator of spatial diffeomorphisms is written as = PmM' + Pfj^R' + PqQ', 
regardless of the particular form of the canonical coordinate transformation. It is thus equating these two expressions 
of the super-momentum Hr, with M, Pm, and Q written as functions of A, R, F and their respective momenta, that 
gives us one equation for the new Pj^ and Pq. This means that we have two unknowns, Pj^ and Pq, for one equation 
only. This suggests that we should make the coefficients of R' and Pp equal to zero independently. This results in 

Pr ^ Pr- lid - 3) fc R-^p-^AP^ - i(d - l)r2pp-iAPA + \{d- 3)if^-2p2^-2rf+5^-i^p^ 
-^R-^APa + R"A-^F-^Pa - A'R'F-^A-^Pa - P^F'^A-^ 

+ {d~3)R'R-^F-^A-^PA, (67) 
Pq = -KdT - K^^R-^'^'^-^^PrP-'APA. (68) 
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(Note that in addition one finds Pq = —KdC-) It can then be shown that Pq commutes with the new Pj^ and with 
the rest of the new coordinates, except with Q. We have now all the canonical variables of the new set determined. 
For completeness and future use, we write the inverse transformation for A and Pa, 



A 



PliF) 



Pa = 2a{d-2)R'^-^FPM {{R'fF-^ - 
In summary, the full set of canonical transformations are the following, 
R = P, 

a{d-2)R'^-^ 



PliF) 



(69) 
(70) 



M 

Q 



RHd-3) 



= PvK- 



Pr = 



PR-^{d-3)kR-'F-^APA 



-(d-l)r2pp-iAPA 



3)K^^P^R-^'^+'^F-^APa 



-R-^APa 



A'R'F-^A-^Pa - P'.F-^A-^ 



-A + P"A-ip-iPA 

+ {d-3)R'R-^F-^A-^PA, 
i2a{d-2))-^R-^'^-^^F~^APA, 

-KdT ~ K-^R-^^'^-^^PpF-'^APa. (71) 

It remains to be seen that this set of transformations is in fact canonical. In order to prove that the set of equalities 
in expression l(7T|) is canonical we start with the equality 

2a{d-2)R'^-^R' + APa 



Pm 

Pq 



Pa6A + PrSR + PrSV 



PmSM - P^SR - PqSQ = [a{d- 2)R'^-^6Rln 

2a{d-2)R 



(5 rPr + APa + a{d - 2)R'^-'^R' In 



2a{d - 2)R'i-^R' - APa 
'^-^R' - AP/ 



+ 



2a{d-2)R'i-3R' + APa J ' ^"^^^ 

We now integrate expression (|72|) in r, in the interval from r = to r = co. The first term on the right hand side of 
Eq. ll72]) vanishes due to the falloff conditions (see Eqs. (I25|-(l33| and Eqs. (l34])-(|42l). We then obtain the following 

expression 



/ dr (PaM + PrSR + PrST) dr {PmSM + P^5R + Pq5Q) = 5lo [A, P, T, Pa, Pr] , 

Jo Jo 



(73) 



where 5uj [A, P, F, Pa, Pr] is a well defined functional, which is also an exact form. This equality shows that the 
difference between the Liouville form of {P, A, F; Pji, Pa, Pr} and the Liouville form of {R, M, Q; Ppj^, Pm, Pq} is 
an exact form, which implies that the transformation of variables given by the set of equations iffl]) is canonical. 

Armed with the certainty of the canonicity of the new variables, we can write the asymptotic form of the canonical 
variables and of the metric function F{t, r). These are, for i — > 

ARl{t)A^^(t)r^+0{r^), 
Ro{t) + R2{t) r^ + 0{r^), 
Mo{t) + M2(t) r^ + Oir^), 



Fit,r 
R{t,r 
M{t,r 
Q{t,r 
PR(t,r 
PM{t,r 
PQit,r 



Qo{t)+Q2{t) r''+0{r^)., 

0{r), 

0{r), 

0{r). 



with 



Mo 



a(d 
a{d - 



2)R 



d-3 



-2 n2 



k + Pf XT^P, 



2p-2(d-3) 



(74) 

(75) 
(76) 
(77) 
(78) 
(79) 
(80) 

(81) 



2) (d-3)Pg 



rf-4 



r2p2 



P, 



■2{d-3) 



P2 



+P^' 



2Q0Q2 2{d~3)QlR2 



P, 



2(rf-3) 




P; 



2d-5 



21 R0R2 



4P2A-2 



(82) 
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For r ^ oo, we have 



R{t,r) 


= r + rp(t)r ^ ' 4 


- <J [r 


M{t, r) 


— M^[t) + O (r ^ 




Qit,r) 


= Q+(t) + 0°°(r-('^- 










PM{t, r) 






PQit,r) 







(83) 
(84) 
(85) 
(86) 
(87) 



where M+{t) = 2a{d- 2) (A(i) - {d - l)p{t)), as seen before in Eq. (06]) and Eqs. ([63l) - l[64l) . 

We are now almost ready to write the action with the new canonical variables. It is now necessary to determine 
the new Lagrange multipliers. In order to write the new constraints with the new Lagrange multipliers, we can use 
the identity given by the space derivative of M, 

M' = -K-^ [r'H + (2 a{d - 2))-^R-^'^-^^P!y {Hr - TG)) +2a{d- 2) KjPrPl^R-^'^-^^ . (89) 

Solving for H and making use of the inverse transformations of A and Pa, in Eqs. i(69l) and ((70l) . we get 

M'F-iR' + FPMPR + 2a(d- 2)i?-(''-3)QQ/^/^-i 
H = 1 , (90) 

Hr = PmM' + P^R' + PqQ' , (91) 
G = -KdQ' . (92) 

Following Kuchaf the new set of constraints, totally equivalent to the old set H{t,r) — 0, H.r(t,r) = 0, and 
G = 0, outside the horizon points, is M'{t,r) = 0, Pji{t,r) = 0, and Q'{t,r) = 0. By continuity, this also applies on 
the horizon, where F{t, r) = 0. So we can say that the equivalence is valid everywhere. 
The new Hamiltonian, which is the total sum of the constraints, can now be written as 

NH + N^'Hr + $G = N'^'M' + N^P^ + N^Q' . (93) 

In order to determine the new Lagrange multipliers, one has to write the left hand side of the previous equation, Eq. 
((93l) . and replace the constraints on that side by their expressions as functions of the new canonical coordinates, spelt 
out in Eqs. l(90l) - |92|) . After manipulation, one gets 



NM ^ -NP-^R'A-^ + {2aid-2))-^N''R-'^'^-^^F-^APA, (94) 
^ -{2a{d~2))-^NR-'^'^-^^Pf, + N'^R' , (95) 
= 2a{d-2)NR-'^'^-^'>A-^K^^PrF-^R' - N^'KdT - N^'K^^R-^^'^-^^PrP-^APA, (96) 

allowing us determine its asymptotic conditions from the original conditions given above. These transformations are 
non-singular for r > 0. As before, for r ^ 0, 

N''{t,r) = -^N^{t)AoR^^ + 0{r^) , (97) 

N^{t,r) = 0{r^), (98) 
iV'3(i,r) = ~KdMt)+a{d-2)NiQoAoR^^Ro'^''-^^ +Oir^), (99) 

N^{t,r) = -iV+(i)+0°°(r-('^+i)), (100) 
7V^(t,r) = 0°=(r-('^-2))^ (jQj) 
N^{t,r) = -Kd^+{t) + 0°°{r-^'^-^'>). (102) 

The conditions l(97|) - l|102p show that the transformations in Eqs. (|94l) - |95l) are satisfactory in the case of r ^ oo, but 
not for r — > 0. This is due to fact that in order to fix the Lagrange multipliers for r ^ oo, as we are free to do. 



and for r ^ oo we have 
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we fix N+{t), which we already do when adding the surface term — / dtN+M+ to the action, in order to obtain the 
equations of motion in the bulk, without surface terms. The same is true for $+. However, at r = 0, we see that 
fixing the multiplier 7V*^ to values independent of the canonical variables is not the same as fixing Nik^^ to values 
independent of the canonical variables. The same is true of the fixation of N'-^ with respect to $o- We need to rewrite 
the multipliers and N'^ for the asymptotic regime r ^ without affecting their behavior for r — > oo. In order to 
proceed we have to make one assumption, which is that the expression given in asymptotic condition of AI{t,r), as 
r ^ 0, for the term of order zero Mq{Rq, Qq), defines Rq as a function of Mq and Qo, and Rq is the horizon radius 
function, Rq = R^^Mq^Qq). Also, we assume that Mq > M^^h{Qq), where M„ij((5o) is found through the system 
of two equations R'^^'^'^^ F{R) = and {R'^^'^'^^ F{R)y 0, with F{R) defined in ©. With these assumptions, we 
are working in the domain of the classical solutions. We can immediately obtain that the variation of Rq is given in 
relation to the variations of Mq and Qq as 



SRq = {a(d - 2) ((d ~3)k R^-^ + {d- l)r^R^-^ - (d - 3)Q^i?-'''~')) } ' 



- ^,Mo - ^-(1;J)Q° ^Q„ ) . (103) 



r: 



This expression will be used when we derive the equations of motion from the new action. We now define the new 
multipHers 7V*^ and N'^ as 



= -N 



M 



(1 - <?) + 2<-3^ ((d -3)k + id- l)r _ _ 3)q2^- 



(d-2) 



(104) 
(105) 



where g{r) = 1 + O(r^) for r ^ and g{r) = 0°°{r ^) for r — > oo. The new multipliers, functions of the old 
multipliers 7V*^ and N'^ , have as their properties for r ^ 0, 



and as their properties for 



00, 



N^'\t,r) = iV+(t) + 0°°(r-(''+i)), 



(106) 
(107) 



(108) 
(109) 



When the constraints M' = — Q' hold, Nq is given by 



M 



(110) 

With this new constraint 7V^'^, fixing A^iA^^^ at r = or fixing Nq'- is equivalent, there being no problems with n'^, 
which is left as determined in Eq. (|95l) . With respect to iV"^ the same happens, i. e., fixing the zero order term of the 
expansion of iV*3 for r ^ is the same as fixing $o, even if multiplied by a dimension dependent constant K^. At 
infinity there were no initial problems with the definitions of both N^^ and N'^ . 

The new action is now written as the sum of Ss, the bulk action, and Sos, the surface action, 

s[m,R,Q,Pm,Pji,Pq;N^',n'^,nQ = 

/nOO 
dt J dr |PmM + Pj^R + PqQ + N'^Q' - N^P^^ + N^^{1 - g) M' 

+iV*^ 2 g ((d - 3) fc R^-^ + {d- 1)1-^ R^-^ ~ {d - 3)QgP^<'^"'^) 
X {R^-^M' -2a{d-2)Q'Qo)} 

dt{(2a N^' R^-^ - N+M+^ + {^oQo - } ■ (111) 

The new equations of motion are now 

M = 0, (112) 
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R 




(113) 


Q 


= 0, 


(114) 


Pm 




(115) 


A? 


= 0, 


(116) 


Pq 


— ( ArQ\' 


(117) 


M' 


= 0, 


(118) 




= 0, 


(119) 


Q' 


= 0, 


(120) 



where we understood TV*^ to be a function of the new constraint, defined through Eq. I|104p and N'^ as a function of 
the new constraint defined through Eq. (|105p . The resulting boundary terms of the variation of this new action, Eq. 
(|llip . are, first, terms proportional to SM, 6R, and SQ on the initial and final hypersurfaces, and, second. 



dt ( 2 aR^-^6N^'' - M+SN+ ] + Kd { QoS^o - Q+S^+ 



tM 



(121) 



To arrive at l|12ip we have used the expression in Eq. (|103p . The action in Eq. (|llip yields the equations of motion, 
Eqs. I|112p - (|120p . provided that we fix the initial and final values of the new canonical variables and that we also fix 



the values of Nq^ 

'tM 



and of Nj^, and of $o and Thanks to the redefinition of the Lagrange multiplier, from N to 



, the fixation of those quantities, Nq and A''+, has the same meaning it had before the canonical transformations 



and the redefinition of N . The same happens with $o and This keeping of meaning is guaranteed through 
the use of our gauge freedom to choose the multipliers, and at the same time not fixing the boundary variations 
independently of the choice of Lagrange multipHers, which in turn allow us to have a well defined variational principle 
for the action. 



D. Hamiltonian reduction 



We now solve the constraints in order to reduce to the true dynamical degrees of freedom. The equations of motion 
(|112p - l|120p allow us to write M and Q as independent functions of space, r. 



M{t,r) 
Qit,r) 



mit), 
q(t). 



The reduced action, with the constraints taken into account, is then 



S 



m, Pm, q, Pq; K'\ ^o, *+ 



= / dt 



Pmrii + Pqq - h j 



where 



Pm 



Pq 



M ■ 



/•oo 

/ drPQ, 
Jo 



and the reduced Hamiltonian, h, is now written as 



h(m, q; t) = -2a N^' R^^ + iV+m + lU q ( $+ - $o 



(122) 
(123) 



(124) 

(125) 
(126) 

(127) 



with being the horizon radius. We also have that m > Af„,it(q), according to the assumptions made in the 
previous subsection. Thanks to the functions NQ^{t), N+{t), $o(0, and <&+(i) the Hamiltonian h is an explicitly time 
dependent function. The variational principle associated with the reduced action, Eq. (|124p . will fix the values of m 
and q on the initial and final hypersurfaces, or in the spirit of the classical analytical mechanics, the Hamiltonian 
principle fixes the initial and final values of the canonical coordinates. The equations of motion are 



m = 0, 



(128) 
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Pm = 2a{d-2)N^R^-''{a{d-2)[{d-3)kR'^-^ + {d-l)l-^R^-^-{d-3)q'R- 

{a(d - 2) (d - 3) fc + (d - l)r2i?^-2 _ _ 3)q2i?j^ 



Pq = -(2a(d-2))2qiV, 



M 




(d-2) 



} 



-1 



(129) 
(130) 

(131) 



The equation of motion for m, Eq. I|128p . is understood as saying that m is, on a classical solution, equal to a function 
of the mass parameter of the solution, m = 2a{d — 2)M, where the solution is given in Eqs. ([T]) and The same 
goes for the function q, where Eq. (|129p impHes that q is equal to the charge parameter on a classical solution, Eqs. 
^ and (IH). In order to interpret the other equation of motion, Eq. (|130p . we have to recall that from Eq. (|65l) one 
has Pm — —T', where T is the Killing time. This, together with the definition of Pjn, given in Eq. I|125p . yields 



(132) 



where Tq is the value of the Killing time at the left end of the hypersurface of a certain i, and is the Killing time 
at spatial infinity, the right end of the same hypersurface of i. As the hypersurface evolves in the spacetime of the 
black hole solution, the right hand side of Eq. (|129p is equal to Tq — T^. Finally, after the definition 



Pq = Kd (Co - C+) 



(133) 



obtained from Eqs. (|49l) . (|68|) . and l|126p . Eq. (|13ip gives Pq oc — which is the difference of the time derivatives 
of the electromagnetic gauge ^(t, r) at r = and at infinity. 



E. Quantum theory and partition function 



The next step is to quantize the reduced Hamiltonian theory, by building the time evolution operator quantum 
mechanically and then obtaining a partition function through the analytic continuation of the same operator [ly|. The 
variables m and q are regarded here as configuration variables. These variables satisfy the inequality m > M„,it(q). 
The wave functions will be of the form ■(/'(m, q), with the inner product given by 



(V',x) 



/idrndq V'X j 



(134) 



where A is the domain of integration defined by m > M„,it(q) and /x(m, q) is a smooth and positive weight factor for 
the integration measure. It is assumed that /j, is a slow varying function, otherwise arbitrary. We are thus working in 
the Hilbert space defined as L^{A; ^dmdq). 

The Hamiltonian operator, written as h(i), acts through pointwise multiplication by the function h(m, q;t), which 
on a function of our working Hilbert space reads 



h{t)ip{m, q) = h(m, q; t)i/'(m, q) . 



(135) 



This Hamiltonian operator is an unbounded essentially self-adjoint operator. The corresponding time evolution 
operator in the same Hilbert space, which is unitary due to the fact that the Hamiltonian operator is self-adjoint, is 



K{t2;ti) = exp 



dt' h(t') 



(136) 



This operator acts also by pointwise multiplication in the Hilbert space. We now define 



r := 


r dtN+{t), 


(137) 


e := 


r dtN^'{t), 


(138) 


So := 


/ dt^oit), 


(139) 




r dt^+{t). 
Jti 


(140) 
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Using (|127p . I|136p . (|137p . I|138p . I|139p . and (|140p we write the function K, which is in fact the action of the operator 
in the Hilbert space, as 

K(m;T,e,So,5+) = exp -imT + 2mi?J^~2e-iXdq(S+ - So) . (141) 

This expression indicates that K(t2;^i) depends on ti and t2 only through the functions T, 9, Sq, and S+. Thus, 
the operator corresponding to the function K can now be written as isr(T, 0, Sq, 5+). The composition law in time 
'K(t^;t2)^{t2]ti) — K(t3;ii) can be regarded as a sum of the parameters T, O, Sq, and S-|_ inside the operator 
K(T, 0, So, S+). These parameters are evolutions parameters defined by the boundary conditions, i.e., T is the 
Killing time elapsed at right spatial infinity and 8 is the boost parameter elapsed at the bifurcation circle; Sq and 
S-(_ are line integrals along timelike curves of constant r, and constant angular variables, at r = and at infinity. 

IV. THERMODYNAMICS 
A. Generalities 

We can now build the partition function for this system. The path to follow is to continue the operator to 
imaginary time and take the trace over a complete orthogonal basis. Our classical thermodynamic situation consists 
of a d-dimensional spherical, planar, or hyperbolic charged black hole, asymptotically AdS, in thermal equihbrium 
with a bath of Hawking radiation. Ignoring back reaction from the radiation, the geometry is described by the 
solutions in Eqs. ([IJ-lISl) and |[5])-|[6l). Thus, we consider a thermodynamic ensemble in which the temperature, or more 
appropriately here, the inverse temperature (3 is fixed, as well as the scalar electric potential 4>. This characterizes a 
grand canonical ensemble, and the partition function -Z(/3, 0) arises naturally in such an ensemble. To analytically 
continue the Lorentzian solution we put T = — i/3, and Q = — 27ri, this latter choice based on the regularity of the 
classical Euclidean solution. We also choose So = and S_|_ = il3(j). 

We arrive then at the following expression for the partition function 

Z(/3, (j)) = Tr \k{-i(3, -2m, 0, i^0)j . (142) 

From Eq. I|14ip this is realized as 

Z(/3,0)= j udmdqcxp -/3(m - Xd q0) + 4 a 7ri?J^"2 (in|m) . (143) 

Since (m|m) is equal to (5(0), one has to regularize (|143p . Again, following the Louko- Whiting procedure we 
have to regularize and normalize the operator K beforehand. This leads to 

ZreniP, (fy ^ M j rfmdq exp -/3(m - i^d q0) + 4 a 7ri?j^"^ , (144) 

where A/" is a normalization factor and A is the domain of integration. Provided the weight factor /i is slowly varying 
compared to the exponential in Eq. I|144p , and using the fact that the horizon radius is function of m and q, the 
integral in Eq. (|144p is convergent. Changing integration variables, from m to where 

m = a{d- 2)i?^-3 (r2i?2 + fc + c^K^^'^^'^^) , (145) 

the integral Eq. (|144p becomes 

2ren(/3, cb) ^ M [ JldR^dq exp(-/,) , (146) 
Ja' 

where A' is the new domain of integration after changing variables, and the function /*(i?jj, q), a kind of an effective 
action (see 0), is written as 

d-3fi-2T,2 , , , „2n-2(d-3)\ nr^.^i A „ ^ T)d-2 



h{R^,q):^Pa{d~2)R'^'[l-'R^ + k + ci'R^'''' ~ (3KdCi(j) ~ AaT:R'^\ (147) 

The domain of integration. A' , is defined by the inequalities < and q^ < i?^''* '^^ {k + ^^^~^^h)' '^^^ 
weight factor includes the Jacobian of the change of variables. Since the weight factor is slowly varying, we can 
estimate the integral of -Zren(/3, 4>) by the saddle point approximation. 
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B. Effective action and critical points (i.e., the solutions of the system) 



Now we calculate the critical points of the effective action (|147p in order to evaluate the integral of the partition 
function (|146p through the standard saddle-point method. We need the critical points because the saddle point 
method requires the Taylor expansion of the effective action around one of them. The critical points are found by 
taking the first derivative of the effective action with respect to and to q and making them zero. With this, the 
Taylor expansion in the first three terms has only the zeroth order term, the effective action evaluated at the critical 
points, and the second order terms evaluated at the critical points as well. Higher orders are ignored. This evaluation 
is done in a close neighborhood of a critical point, thus making the zeroth order term more important than the second 
order term. Which one of the critical points is chosen for this approximation is seen below. 

The critical points are found through equating to zero the first derivatives of the effective action l|147p , with respect 
to and to q, i.e., 

= = (148) 



9q 



Two distinct sets of critical points are then found: 

(i) The critical points are given as a pair of values {R^ , q) 



= k Z '}', , (150) 



2a{d-2) 

where, as a reminder, Ka = 4a^27r (d-2)(d-3). 
(ii) Equation (|147p also has a critical point at 



R^^ - 0, (151) 
q = 0, (152) 



which for d > 4 is equivalent to (i?jj,q) = (0,0). For d — A this critical point does not exist, but in practice in the 
study of global minima this makes no difference. 

The critical points belong to the domain of the effective action (|147p , and so the effective action has a derivative 
equal to zero at these points in Eqs. (|149p - l|152p . From (|149p . if 

one finds there is at least one critical point. In more detail we find the following: if in Eq. I|153p there were a < instead 
of a >, then there would be no critical points. If the equality is satisfied in Eq. I|153p . then there is only the critical 
point given by — {2TTp)/{{d — l)/3) and the corresponding value of the charge, by replacing in Eq. I|150p . In 
this last case ~ "^h' If the inequality holds in Eq. I|153p . then there are three situations (a) K^cj)"^ < k{2a{d — 2))^, 
(b) ifJ(/>^ = fc (2 a{d - 2))2, and (c) K'^cj)'^ > k{2 a{d - 2))^. In (a) there are two critical points, given by (i?^, q)^ in 
(|149p - l|150p . in (b) there are two critical points given by R^ = (47rP)/((d— l)/3) and R^ — 0, with the corresponding 
q* given in l|150p . and finally in (c) we have only the critical point given by (i?j^, q+), because R^ < is not physically 
relevant. All this discussion neglects the critical point at (0, 0) which is independent of the value of 0, which must 
be added to the discrimination above. One must also note that the topology was not taken into account in detail. 
If the topology is taken into account, one sees that for k = 0,-1 there is always only one critical point of physical 
significance, given by {R^,<1^) as (jP is positive and, for k = 0,-1, the right hand side of Eq. I|153p is negative. In 

this case the critical point at (i?j^,q~) is outside the domain of physical interest as R^ < 0. 

We can now write the action evaluated at the critical points, where also /3 and (j) are determined as functions of the 
critical pair (i?^,q±). The full expression is 

2) ( kRf'-'^-d^'-l-^Rf-'' 

" \[d-i)k Rf'-'' -id- 3)q±^ + (d - ' ' ^'''^ 
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where a = E^_2/167r, see Eq. (fTOl) . Expression l|154p is quite general. It can be compared with other studies that 
have chosen a particular dimension of spacetime, or a particular asymptotic regime, or a particular topology. 

For instance, for generic d and k — 1 one finds from (|154p an action for d-dimensional Reissner-Nordstrom black 
holes which was not explicitly shown in [l^ . [T3 |. If, for generic d and fc = 1, we further put q = one finds from 

CH the following action h = A.aRf-'\R^''-'-'^-l-^Rf '-'^)/[{d~i)Rf '-'^ + {d-l)l-^^^^^^^ 
which is the action found in In addition, by putting d — 4 and k — 1 into l|154p one finds /* = 

7ri?^ (^R'^ — — / (^h ~ Sl^^R^ which is the action first found in If we further assume q = 0, 

one finds the Hawking-Page action = ttR^ (^R^ ~ l^'^Rf^ I + 3^"^^?^) [8]. Or, by putting d = 4 and A: 

into l|154p one finds /, = 7r_R^ ^q^ + / — "il'^R^ , which was found in |ilO|]. We defer a full comparison 

to the next subsection, after having studied the phase transitions of the system. 

C. The most stable solutions and phase transitions 

1. Analysis in d dimensions 

Now, the ensemble in this case has already been identified as the grand canonical, where T (or p) and (j) are fixed. 
What now is needed is to obtain the most stable thermodynamic solutions for this ensemble. In order to do that, 
it is necessary to know the interval of the order parameter in which a given critical point is a global minimum, 
local minimum or none of the latter. Each different minimum represents a physical system, such as hot fiat space, for 
instance. Note that, hot fiat spaces and black holes, are two different sectors of the solution space, one with trivial 
topology, the other with black hole topology, respectively. To be definitive, let us fix the temperature, or /3, on the 
boundary. Now, (j) is also fixed on the boundary, but we can fix it with any value we like. So, imagining changing the 
parameter 4>, fixed on the boundary, one finds that, for instance, one local minimum may turn into a global minimum, 
signaling that a phase transition occurs at a certain value of 4>. 

The saddle point approximation of the partition function l|146p chooses the global minima of the critical points for 
each value of 0, on the understanding that (3 is fixed. From Eqs. (|149p . (|15Qp . I|15ip . and (|152p one finds the values 
of the critical points in terms of and q. Eqs. I|15ip and l|152p are at the origin of the {R}^, q), and do not depend 
on the parameter 0. Eqs. (|149p and (|150p will change in value as 4> changes. The expression (|150p implies that once 
is determined for a critical point, q is immediately given. So, on the plane (i?jj,q), the critical points are found 
by the curve defined through l|150p . for each given (f). Then, if we replace the value of the charge as a function of 
i. e., q(i?jj) as given in l|150p . into the effective action (|147p we obtain a one variable function of written as 

= Pa{d-2)(l-'^R^^ + kR^^ + Kl(t)'^R^^{2a{d-2))-^^- 
|3{2a{d - 2)y^Kj(j)^R^^ - AnaR^-'^ , (155) 

which is a function of Now, not all R-^ are solutions, i.e., critical points, only those defined by l|149p belong to 
the solution manifold. Choosing from those that belong to the solution manifold we are interested in solutions that 
minimize the effective action. In order to find out which is a minimum, for a given cj), and even which is the lowest 
minimum of them all, it suffices to analyze the one variable function given in (|155p . This function could have been 
written with respect to q, thus allowing us to find the minima in q. However, as we are dealing with a black hole 
whose entropy we also want to determine, and as the entropy is proportional to the horizon area, which is itself a 
function of R^, it is much more convenient to study the effective action as a function of R^. Now, the first critical 
point is at = 0. At this point the action l|155p is equal to zero, /*(0) = 0. This critical point is always there, as 
= does not depend on any one parameter, except perhaps the dimension, which is considered here as d > 4. 
If d = 4 then R^ = is no longer a critical point, but for any other critical point to be a global minimum it will 
have to obey that /*(i?^) < 0. This is valid for all the other dimensions where = is in fact a critical point. If 
/*(i?j^) > 0, the critical point in question is not a global minimum. All this reasoning is valid when there are no two 
critical points such that /*(i?j^) < 0. In the case of two critical points with /,(i?j^) < 0, the one that gives the least 
value of l|155p is the global minimum. So the first step is finding the zeros of (|155p . Not only do the zeros help finding 
the global minimum, but also they give the intervals for the order parameter for each different global minimum. In 
other words, the zeros help discriminate the physical phases. The zeros of the effective action (|155p are given by 



{Rir-' = 0, 



(156) 
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The first zero, = is independent of (/). In relation to the second zero, there follow different conclusions, depending 
on the value of (j)'^ for . Indeed, R*^ must respect 

Kj(f>^ > -16oV/2/3-2 + fc(2a(d-2))2, (158) 

so that at least a zero of the form of (fTSTl exists. Then, if ifj^^ were smaller than the right hand side of (|158p there 
would be no new zeros, besides R'^. This would mean physically that there would be no black hole, and there would 
only be a heat bath with temperature (see [8] for the particular case of Schwarzschild-AdS). If the equality of l|158p 
holds, then we have a double zero at — {2nP)/{{d— 2)(3). This point is also a critical point, i.e., it is a solution 
for the system. If we replace the equality in l|158p back in l|149p . we see that R^ = {2nP)/{{d — 2)(3). Here both 
= and — {2iT P) / {{d ~ 2)j3) are global minima, so there is a degeneracy. Here, there is an equal probability 
for the system to consist of a heat bath alone or a charged black hole of radius = {2iTp)l{{d— 2)(3) with thermal 
radiation. The equality in l|158p marks the threshold from which the critical point is a global minimum. Once the 
inequality in l|158p holds, then is the only global minimum, ending the degeneracy, i.e., I^,{R^) < /*(0) = 0. This 
implies the stable solution consists of a charged black hole of radius R^ in equilibrium with a heat bath temperature 
f3~^. Perhaps some other minor comments in relation to R^ are in order. The other critical point R^ is always 
less than R^, and the effective action l|155p is always positive at R^ as long as i^J^^ — k{2a{d— 2))^ > 0. Once 
if^i^^ — fc (2 a(d—2))2 = 0, then both the critical point = and /*(i?j^) = 0. At the same time i?^ =0. However, 
the action l|155p at is always larger than at R^, once there are two different critical points R^. If any critical 
point R^ is less than zero, it has lost its physical meaningfulness. 

All the previous discussion is general, and applies directly to the spherical topology. However, if the topology is 
different, i.e., k = 0,-1, then the inequality in the expression l|158p holds always. This means that R^ is always the 

global minimum. Therefore the classical solution is found in the global minimum of the critical point (i?J^, q"*") of the 
effective action /*(i?jj, q) (Eq. (|147p ). Differently from the spherical case, there are no phase transitions, as discussed 
in 

For an illustration of the previous discussion, a choice is made, namely, d = 5, A: = 1, a = 1, and I = -\/T0. We also 
choose (i = ^tt. Then, the effective action is written in the form 

Tr~^h{RY,) = - 4i?3 + ci?^ , (159) 

where c = 10(1 — ^ncfP'). So, for a choice of c = 0, 3, 4, 4.4, 4.5, 5 respectively, which amounts to a choice of different 
values of 4>, the order parameter for this phase transition, we plot, in Figure[2l 7r~^/*(i?jj) as a function of These 
plots illustrate the process of phase transition from a phase where the black hole is a stable solution to a phase where 
the action has a global minimum for hot fiat space. In this instance the phase transition happens due to a change in 
the order parameter 0, the electric potential. The evolution in the plots is toward a higher value of c, or smaller value 
of (j). In more detail, in the first plot there is a critical point at a finite radius where the effective action is negative. 
The other critical points are at zero radius. Thus the global minimum, the stable solution, is a black hole of radius 
given by the critical point. This is expected on physical grounds, since a higher 4> means a higher electrical pressure 
of the walls at infinity on the charged particles which then tend to concentrate at the center, forming a black hole. 
The next plot shows two distinct nonzero critical points, plus the critical point at the origin. The one with larger 
is the stable solution, where again the effective action is negative. In the other critical points, the action is either zero 
or positive. The next plot shows us the phase transition. Here the value of c is such that the effective action at the 
larger critical point is equal to zero, as it is zero at zero radius. This implies that the effective action has two global 
minima. Here the black hole with the larger radius and hot fiat space are equally probable. The system is then a 
mixture of two states which can transit from one to the other. The value oi <f)\s (j) = . For higher values of c, the 
likeHest outcome is now hot fiat space, because the effective action is positive for any value of c, except at the origin. 
Even so, the larger radius critical point is still a local minimum, which allows for a metastable black hole solution. 
However, there is a limit for the existence of such metastable solutions. The last but one plot, for c = 4.5, shows the 
merging of the two nonzero critical points into one, where the effective action is positive and where it has neither a 
maximum nor a minimum. It is thus a completely unstable solution. For values of c such that c > 4.5 there are no 
other critical points, except for the zero radius, which is equivalent to saying that hot fiat space is the only outcome 
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FIG. 2: Plots of the effective action, tt ^7*(7?Jj) in five dimensions, for the choices d — 5, k = 1, a = 1, I — \/T0, and 
c = 0, 3, 4, 4.4, 4.5, 5, respectively, where c = 10(1 — ^■K(j)'^), see text for details. 



of the ensemble. There is thus a value of (j> below which there are no black holes for the ensemble. In our example the 

value is = (gff;;:) ^ ■ This is analogous to the Bose-Einstein condensation, which is a phase transition. In this case 
one has as variables the chemical potential /i and the temperature T, or (3. If one fixes T and increases the number of 
particles and so increases /i one has a phase transition to a condensate. So /i is equivalent to here, the difference is 
that bosons feel an effective attractive force, whereas charged particles with the same sign feel a repulsive force. So, 
when one lowers one has hot flat space, whereas when one raises fi one has a condensate. When one raises one 
has a black hole, whereas when one lowers one has a free gas. 

Now, we can instead study the phase transition by fixing (f) a priori, and keeping T , or more precisely, the inverse 
temperature /3, fixed, but allowing it to change from situation to situation. For the same d = 5, fc = 1, a = 1, and 

I = VTO, but now with the value of (j) given a priori by cj) = j\/ , the effective actions reduces to 



7r-i/,(i?i,) = 



3/3 ^4 



27/3 2 
10^ ■ 



(160) 



Through the choice of different values of we can make a series of plots which show the different phases of the 

ensemble with respect to the inverse temperature (3. The choices for are = 11)' "RJ' §' § + 0-01, § + 0.1. 
We plot, in Figure O the effective action 7r~^/*(i?jj) as a function of iJj^, for those values of /3. These plots illustrate 
the process of phase transition from a phase where the black hole is a stable solution to a phase where the action 
has a global minimum for hot fiat space. In this instance the phase transition happens due to a change in the order 
parameter /?, the inverse temperature. The evolution in the plots is toward a higher value of f3, or smaller value of 
the temperature. In more detail, in the first plot there is a critical point at a finite radius where the effective action is 
negative. The other critical points are at zero radius. Thus the global minimum, the stable solution, is a black hole 
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FIG. 3: Plots of the effective action, it in five dimensions, for the choices d = 5, fc = 1, o = 1, Z = VW, and 

^ = ^, ^, |, I + 0.01, I +0.1, respectively, see text for details. 

of a radius given by the critical point. Note that the choice of /3 = would have yielded an effective action given by 
= — 47ri?j^, which would not have given any critical point except = 0. This infinite temperature ensemble 
would have had a global niaxiniuni at zero radius for the effective action, implying that hot flat space at infinite 
temperature is unstable, and at the same time showing that only a black hole of infinite radius would satisfy the 
conditions for stability. The next plot shows two distinct nonzero critical points, plus the critical point at the origin. 
The larger is the stable solution, where again the effective action is negative. In the other critical points, the action 
is either zero or positive. The next plot shows us the transition. Here the value of (3 is such that the effective action 
at the larger critical point is equal to zero, as it is zero at zero radius. This impHes that the effective action has two 
global minima. Here the black hole with the larger radius and hot flat space are equally probable. The system is then 
a mixture of two states which can transit from one to the other. The value of /3 is = 22zL, por higher values of (3, the 
likeliest outcome is now hot fiat space, because the effective action is positive for any value of /?, except at the origin. 
Even so, the larger radius critical point is still a local minimum, which allows for a metastable black hole solution. 
However, there is a limit for the existence of such metastable solutions. The last but one plot, for (3 = iSivI^ gi^Q,j^g 
the merging of the two nonzero critical points into one, where the effective action is positive and where it has neither 
a maximum nor a minimum. It is thus a completely unstable solution. For values of /? such that P > there 
are no other critical points, except for the zero radius, which is equivalent to saying that hot fiat space is the only 
outcome of the ensemble. Again, we can compare this phase transition with the Bose-Einstein condensation. Now, 
one fixes ^ in one case and fj, in the other. So when one raises /? (lowers the temperature) one has hot fiat space, and 
a condensate in the other case. When one lowers /3 (raises the temperature) one has a black hole, and a free gas in 
the other case. 

We note that there is nothing special happening when we take the number of dimensions to be very large, d — > oo. 
All the quantities suffer a change in scale, certainly, but qualitatively the behavior is analogous. 
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2. Comparison with other results 



We can compare our results for d-dimensional, Reissner-Nordstrom-AdS black holes, with either spherical, planar, 
or hyperbolic horizon topology, with results from other authors in several other instances. We divide this comparison 
into the works that also used Hamiltonian methods and the works that used Euclidean path integral methods. As we 
have seen, Hamiltonian methods perform a Hamiltonian reduction in the Lorentz theory, and then take the trace of 
an analytically continued evolution operator. Euclidean path integral methods perform a Hamiltonian reduction in 
the analytically continued action. The boundary conditions in the two approaches are identical, and the difference 
between them is in the order of quantization (first in the Hamiltonian methods) and Euclideanization (first in path 
integral methods). 

(a) Hamiltonian thermodynamic methods 

Firstly, Reissner-Nordstrom-AdS black holes in four dimensions (d = 4) with spherically symmetric (k — 1), were 
studied by Louko and Winters-Hilt in [131, using a Hamiltonian thermodynamics formalism. They analyzed these 
systems in both the canonical ensemble and grand canonical ensemble. Now, if we put d = 4 and fc = 1 in our action 
action (|147p and in the critical point (|149p - l|150p . it reduces then to the grand canonical ensemble effective action and 
the critical point of [23|. There is a small difference in the coefficient Kd between our work and when d — i, due 
to a different choice in the coefficient of the Maxwell term in the action ^ . Here we have followed Myers and Perry 
[30| . where the Maxwell term is jFf^^F'^'^, whereas in [20] there is no j term. 

Secondly, Reissner-Nordstrom-AdS black holes in four dimensions {d = 4) with planar and hyperbolic symmetry 
(fc = 0,-1), were studied by Brill, Louko and Peldan in f2l| in both, the canonical ensemble and grand canonical 
ensemble, using implicitly a Hamiltonian thermodynamics formalism. When we put d — A and (fc = 0,-1) in our 
action l|147p and critical point l|149p - (|150p one recovers the results in [2l| for the grand canonical ensemble. Here, as 
with the case above, there is the caveat of the constant Kd, i.e., it will not render the same numerical coefficient due 
to the choice of a different constant for the Maxwell term in the action ^ . 

Thirdly, it is interesting to note that the calculations done, also using a Hamiltonian thermodynamics formalism, 
by Louko, Simon, and Winters-Hilt in [l^l and ours match in the appropriate limits. In [23| it was studied a spherical 
black hole (fc = 1) in a finite box with radius tq in a Gauss-Bonnet (Lovelock) theory with parameter A, without 
cosmological constant, in five dimensions, d = 5. Here we have studied charged black holes in general relativity with 
cosmological constant in d dimensions. Thus, in order to match both results, we do the following: in [2§| we take the 
Gauss-Bonnet parameter equal to zero (A = 0) to recover general relativity, and take the box to infinity (re oo); 
here we take d = 5, the cosmological constant and the electrical charge equal to zero (Z — > cx) and q = 0), and assume 
spherical topology (fc = 1). Then one gets as the equation for the critical point of the effective action /3 = 27ri?j^, both 
in jl^l (see Eq. (5.3) of [22]), and here (take the appropriate limit in l|149p ). 

(b) Euclidean path integral methods 

Firstly, the thermodynamics and phase transitions of the spherical Reissner-Nordstrom black holes in AdS space- 
times in d-dimensions can be studied if we put fc = 1 in Eq. I|154p . yielding 

h \ (d - 3) i?f ^'-'^ -{d- 3)q±^ + (d - 1)1^^^^"-^^ j 

the case studied in (l^ . [l3| but in which the action (|16ip was not explicitly displayed. If we further put q = in Eq. 
(|154p . we obtain the action for Schwarzschild black holes in AdS spacetimes in d-dimensions, i.e., 

^ V (d - 3) ^ + {d- ^ ' ) 

This action was found by Witten [ll| through mixed methods. On further choosing d = 4 one finds the Hawking-Page 
action 0. 

Secondly, Reissner-Nordstrom-AdS black holes in four dimensions (d — 4) with spherically symmetric (fc = 1), were 
studied by Pega and Lemos in P| using instead Euclidean path integral methods in the grand canonical ensemble, 
and where the black holes are put inside a rigid box of radius re. As already remarked in the results match, as 
they should, those of jl^l, when tb oo. It is interesting to further compare [§| with our results here. In [§], when 
taking the Hmit re ^ oo, it was found there are two sets of solutions, which are obtained through different ways: (i) 
fix r_|_, which is our R-^, and the charge of the black hole e, which is our q, or (ii) fix the boundary conditions, which 
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means fixing j3 and 4>. This second set yields, as tb ^ oo, a divergent ^ tb and a divergent q ~ r^, and also 
leads to an entropy going to infinity as r|. In this limit the heat capacity and the thermal energy diverge too. These 
very interesting solutions, first found for the asymptotically fiat Schwarzschild black holes, in represent metastable 
black holes which can quantum tunnel into hot fiat space and vice versa. Since we have taken from the start that our 
box is AdS infinity (i.e., tb ^ oo at the outset), our boundary conditions do not recover these metastable solutions, 
and we will not discuss them anymore (see for further analysis of these solutions). The first set, which fixes 
and q, is of our concern now. In f§|, the temperature T and electrostatic potential tend to zero as r^^, and the 
energy tends to zero as well, E = M/ (1 + r^/P)^ ■ Here we have fixed precisely the temperature at infinity (following 
[20| ') so that E = M . But the difference between the two adjustments is merely one of choosing the zero point. After 
these remarks we can proceed, and by putting d 4, fc = 1 into (|154p one finds 



_ ^Rl{Rl-^^-l-^Ri 



Rl - q2 + 3Z-2i?4 



(163) 



"-h ^ ^^h 

which is precisely equal to the one first found in If we further put q = 0, one finds again the Hawking-Page 
action @| /* = 7ri?? + 31 ^Ry^i an action which has attracted much attention. All these black hole 

solutions in AdS spacetime are stable. We can also find the action in the asymptotically fiat case, i.e., I ^ oo. Then, 
for d = 4, fc = 1, and ? ^ oo in Eq. I|154p yields /* = ""^^1 which were first found by York [3| for the uncharged case, 
and by for the charged case. One also finds that these asymptotically fiat black holes are unstable in the grand 
canonical ensemble, although the branch of large black holes not recovered by our boundary conditions are stable (see 
again [1,01, see also for the discussion of the metastability of these objects). 

Thirdly, Reissner-Nordstrom-AdS black holes in four dimensions {d = 4) with toroidal planar topology {k — 0), 
were studied by Pega and Lemos in using Euclidean path integral methods in the grand canonical ensemble, and 
where the black holes are put inside a rigid box of radius vb- As already remarked in the results match, as they 
should, those of fip . for all tb, including when tb — > oo. It is interesting to further compare with our results 
here. Again, in [10||, when taking the limit tb ^ oo, it was found there are two sets of solutions, which are obtained 
through different ways: (i) fix R-^ and the charge of the black hole q, or (ii) fix the boundary conditions, which means 
fixing P and (p. This second set yields, as tb ^ oo, a divergent R-^ ^ tb and a divergent q ~ ^'b, and also leads to an 
entropy going to infinity as r|. In this Hmit the heat capacity and the thermal energy diverge too. As in the spherical 
case, these solutions, represent metastable black holes which can quantum tunnel into hot fiat space and vice versa. 
Since we have taken from the start that our box is AdS infinity (i.e., tb — > oo at the outset), our boundary conditions 
do not recover these metastable solutions, and we will not discuss them anymore (see fl^ for further analysis). The 
first set, which fixes and q, is of our concern now. In the temperature T and electrostatic potential tend 
to zero as r^^, and the energy tends to zero as well, E ^ M/{'KrB/l). Here we have fixed precisely the temperature 
at infinity (following jl^l) so that E = M. But the difference between the two adjustments is merely one of choosing 
the zero point. After these remarks we can proceed, and by putting d = A, k = into (|154p we find 

which was first found in [lo| . 

The analysis displayed in (a) and (b) above, and the comparisons made, show that the new action, Eq. I|147p . for 
general d and arbitrary topologies recovers previous results in the appropriate limits. 



D. Thermodynamic quantities for the solution containing a black hole 



The partition function Zxen{P,<t>) can now be written, for K'^cf)^ > —IQa'^iT^l'^fi ^ + k{2a{d— 2))^, through the 
saddle-point approximation 



lnZren(/3,0) 



P exp 



lnP-/4i?+,q^ 



(165) 
(166) 



where P is a slowly varying prefactor, and /*(i?j^, q+) is the effective action evaluated at the global minimum critical 
point. By ignoring the prefactor's logarithm, which closer to (iJjl^jq"'") is less relevant, we are able to determine the 
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value of m at the critical point, where we find that it corresponds to the value of the mass of the classical solution of 
the black hole given in Eqs. l(I])-ll5]). Thus, when the critical point dominates the partition function, we have that the 
mean energy (E) is given by 



(E) =-—lnZren~a{d-2)R+ " (^r^i?+- + fc + if^ 0^(2 a(d - 2))-^ j = m+ , (167) 
where m+ is obtained from Eq. (|145p evaluated at (^j!j,q^)- The thermal expectation value of the charge is 

(Q) -/3-'^ln2ren«if<iq+, (168) 

as given in Eq. I|150p . We also can write the temperature T = /3~^ as a function of the critical point and the 
electric charge q+ 

id-l)l-^R+' + {d-3)k-{d-3)iq+fR+-"^''~'^ 

T = a T a , 169 

with explicitly written in (|149p . and where q"*" is given in in Eq. I|150p . If we replace the value of the charge 
of the extreme black hole solution, which is given in q^ — R^'^~^^ (jt + , the maximum the charge can 

attain inside the domain of integration of the partition function, into the expression for the temperature, Eq. (|169p . 
then we get a null temperature. It can be shown that dm^/df] < 0, which through the constant (j) heat capacity 
Ccj, — -~0^{d{E) /dp) tells us that the system is thermodynamically stable. 
Finally there is the entropy of the black hole, S, which is given by 

We see that the extreme black hole, despite having null temperature, has an entropy proportional to its horizon area. 



V. CONCLUSIONS 



From a Hamiltonian thermodynamics formalism we have found a quite general effective action l|154p valid for charged 
higher dimensional c? > 4 AdS (negative cosmological constant) spacetimes, with spherical, planar, and hyperbolic 
topologies. There are black hole as well as hot fiat spacetime solutions. It was shown that the phase transitions, 
already present for the Reissner-Nordstrom AdS black hole in four dimensions, show up again for the generic d- 
dimensional charged case. In these phase transitions a black hole turns into hot fiat space and vice-versa. However, 
this phase transition does not happen for nonspherical topologies. We have also calculated the entropy of these 
higher dimensional charged AdS black holes. This entropy is given by a quarter of the horizon area for any topology 
considered. It is thus obtained that the entropy relation with the quarter of the horizon area is kept for arbitrary 
higher dimensions, and that the topology does not affect the relation. Since the Hamiltonian formaHsm yields a correct 
consistent expression for the entropy of a black hole, we have proved that the formalism has a wide applicability, both 
to higher dimensions and to nonspherical topologies. These results may be of value for investigations on AdS/CFT 
physics. 
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